IL NUOVO CIMENTO 



Vol. ?, N. ? 



? 



E. C.G. Stueckelberg: a forerunner of modern physics 

F. ClANFRANl|J%*) AND O.M. LECIAN(tH) (** ) 

PJ) ICRA-International Center for Relativistic Astrophysics 

Dipartimento di Fisica (G9), Universita di Roma, "La Sapienza" , Piazzale Aldo Moro 5, 
00185 Rome, Italy. 



Summary. — The pioneering work of E.C.G. Stueckelberg is briefly analyzed: 
the formalism of the Stueckelberg field, able to describe a massive vector field, 
is reviewed, and some applications are presented. In particular, starting from 
his very first application, devoted to describe nuclear phenomena known in 
the Thirties, later attempts to demonstrate the renormalizability of the model 
are considered. Finally, also string theory and LQG are illustrated to be a 
suitable scenario for the Stueckelberg field. 
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1. Biographical notes 

Ernst Carl Gerlach Stueckelberg was born in Basel on February 1st, 1905. In 1926, he 
got his Ph.D. at Munich under the supervision of Arnold Sommerfeld; then, he qualified 
as a university lecturer at the University of Zurich, till he became Assistant Professor at 
Princeton University. 

In 1934 he provided the first covariant perturbation theory for quantum fields. To quote 
a paper of Lacki et al.pQ, 

The approach proposed by Stueckelberg was far more powerful, but was not 
adopted by others at the time. 

Then in 1935, before Yukawa(2j and by a rather different approach, he proposed to 
explain nuclear interactions as due to the exchange of vector bosons. 
Moreover, the evolution parameter theory he presented in 1941 and 1942 is the basis 
for recent work in Relativistic dynamics. But his great achievement in 1942 was the 
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interpretation of the positron as a negative energy electron traveling backward in time. 
Feynmamnn quoted this result in one of his classic papers [3] . 
Stueckelberg died in 1984 in Basel. 

2. The Stueckelberg field 

StueckelbergJH [5] developed the only up to known formulation of a renormalizable 
theory for a massive Abelian boson (for a recent review see [5]). 

The first model for massive vector particles was the Proca one [8] , which simply produces 
the extension of the electro-dynamics by the introduction of a mass term; in fact, the 
Proca Lagrangian density reads 

(1) CProca = AflF^ + m 2 V^. 

It is clear that such a term provides a violation of the Abelian gauge symmetry, and, after 
several decades, Stueckelberg's work has been recognized as responsible of the renormal- 
izability. After the canonical quantization, one obtains the commutation relations 

(2) 

[V^x);V u (y)} = [V^x);Vj(y)} = 0; [V^x); Vj(y)\ = -iU^ + ^d^^x-y) 
where the function A m satisfies 

(3) {d 2 + m 2 )A m (x - y) = 

After 1945, it became clear that the term ^2<9 M <9„ in the commutation relation ([2]) 
gives rise to ultra-violet divergences, which cannot be eliminated even by the renor- 
malization procedure. However, before the development of the renormalization theory, 
Stueckelberg [4j provided a divergence-free formulation. In his model, the starting point 
is a Fermi- like Lagrangian density for a complex vector field A^, i.e. 

(4) C A = -d»Ald»A v + m 2 A\A»- 
since the Hamiltonian density 

(5) Ha = -d^Ald^A" - m 2 4# 

is not positive definite, one has to impose the analogue of the Gupta-Bleuler condition 
in electro-dynamics, i.e. that the expectation value on physical states of d^A^ vanish 

(6) < phys'ld^A^phys >= 

(a sufficient condition is d^A^ \phys >= 0, being A^ just the positive frequency part). 
In the Proca case, this condition comes directly from equations of motion. 
Unlike QED, relation ([6]) cannot stand, since, from canonical commutation relations, 
which read 



(7) [A fl (xy,A u (y)} = [Al(x);Ai(y)}=0 [A fl (x);Ai(y)}=-ir ] ^A m (x-y), 
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one obtains 

(8) ld»A„(x); d v At(y)} = id 2 A m (x - y) = imA m (x -y)^0. 

Stueckelberg solved this puzzle by introducing a scalar field B(x), whose Lagrangian 
density reads 

(9) L B = d a B^d»B - m 2 B^B, 
with canonical commutation relations 

(10) [B(x);B(y)} = [B^x);B^y)} = [B(x);B^(y)] = iA m (x - y). 

Hence, the consistency condition on physical states, such that the Hamiltonian density 
is positive definite, reads as 

(11) S(x)\phys >= {d^{x) +mB{x)) { - ) \phys > = 

and one can easily demonstrate no contradiction exists with the commutation relations 
0, (fTU|) . Therefore, the full Stueckelberg Lagrangian density is 

(12) Cstueck = -d a A\d»A v + m 2 A\A» + d^d^B - m 2 B^B 
which can be cast in the form 

(13) Cstueck = £proca(W») - (d^ + mB^d^ + mB); 

being W M = A^ — -d u B, it coincides with the Proca Lagrangian density on physical 
states. However, there is a main difference between the two formulations: while the mass 
term in Cp roca destroys the gauge symmetry, that in Cstueck is invariant under Pauli 
transformations, i.e. 

(14) {%Z A B^1 K (di + m 2 )A = 0. 

A kind of invariance is expected to compensate the introduction of the additional field 
B and to lower the number of local degrees of freedom to three. In a physical point of 
view, we can think of the field B(x) as eliminating the scalar term d^A^ of the vector 
field. 

3. The Stueckelberg field and the carriers of nuclear interactions 

The aim of the paper by E.C.G. Stueckelberg was to describe electromagnetic and "nu- 
clear" forces (what we would call electromagnetic, weak and strong interactions) within a 
generalization of the formalism developed for charged particles [7] . This formalism deals 
with the scalar massive field A, which obeys, in presence of matter, the field equation 



(15) 



(<9 M <9 M - I 2 ) A = -4nJ, 
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and will be shown to be equivalent to the retarded-potential method, but will offer the 
advantages of approaching the problem form "gauge" point of view ahead of its tim af 1 )! . 
The Lagrangian density C reads 



(16) 



£ = 



'8tt 



dA dA* 
dx ' dx 



+ 1 A* A 



+ 



A*J + 



dA* 
dx 



S + c.c. 



where A and A* are treated like independent quantities. The "effective" current J e // 
rewrites J e f / = J— ( , s) , as a function of the polarization vector S. From the conjugate 
momentum P — the Hamiltonian density H is found, and, for later purposes, it will 
be expressed as 

(17) H = j d 3 x 



-C + AP + A*P' 



! 



d s x (W + V) = W + V, 



where W = J- 



( OA dA' 
\ dx 5 dx 



I 2 A* A 



+8nc 2 P*P, and V 



A* J 



4ttc(PSo + cc), respectively. Motion equations follow from the introduction of the oper- 
ator K, that allows one to get a straightforward definition of J and S: 



(18) 



h I 1 ' 1 SA* 2 



ox 



(19) (^\ [KjA] = JJL = -^ cSo . 

J and S are functions of the canonical variables p and q, which describe the matter distri- 
bution, and obey the (classical) equations of motion p = {i/K)[K, p] and q — (i/h)[K, q]. 
The quantum theory can be implemented by solving the Schroedinger equation H^(t) = 

ih d ^ d ^ . To this end, the functional ^'(T,t) for the wave function is introduced, such 
that ^'(t,t) ee and, accordingly, the functional K(T,t), such that K(t,t) = K (K 

does not depend on t explicitly). Since the functional must satisfy simultaneously the 
two Schroedinger equations 

(20) C T V'(T,t)= (w + ±-?-)*'{T,t)=0, 



(21) C t *'(T,t) = (V(T,t) + *'(T,t) = 0, 

the wave function ^ is defined by H = W + K. The request that the two Schroedinger 
equations be simultaneously solvable leads to the vanishing commutation relation be- 
tween the operators defined in ([20]) and (f2Tj) . i.e. [Ct, C t ] = 0, from which the expression 
for K 

(22) K'(T, t) = e imt-T)/h Ke -iw(t-T)/h 



( 1 ) Throughout this section,we will maintain the original notation adopted by Stueckelberg, in 
order to appreciate the development of his pioneering intuitions. In particular, the introduction 
of the Stueckelberg held will be understood from a historical point of viewed, i.e. via the 
Dirac-Fock-Podolski approximation, rather than from a modern perspective, as reviewed in the 
previous section. 
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is found; consequently, \l/ admits the formal solution 

(23) #(T,t) = e- lWT/h ^(t), 
where ip(t) satisfies the Schroedinger equation 

(24) K "(T)m = ih^, 

with K"(t) = e lWt / h K e - lWt / h . if matter distribution is described in the configuration 
space, with coordinates {q s }, rather than by means of matter fields, the Hamiltonian 
operator K can be rewritten as the sum of two terms, 

(25) K = J2 K ° = I2( R ° + V ^> 

s s 

where the former depends on the {<? s }'s only, R s = R s (q s ), while the latter is a function 
of both the coordinates {q s } and the field A, V s = V s (q s , A((x),t)) = V s (t), as it will be 
explained in the following. Eq. (|24|) now reads 

<*> (**+tI)*-( j,+ " + tI)*-"- 

where the sum (125[) is taken into account. This Schroedinger equation contains the 
"current term" V s , which contains, on its turn, the field A: assuming that this term is 
proportional to a small number, a series expansion will be performed in order to obtain 
the approximated expression for ip, i.e. ip = ip + ip 1 + ip 2 + 

(27) (VtItV = 0, 



dt 



(28) t R+ !i®\^ +Vi ,° =0 



i dt / 

and so on. Collecting the terms for the proper approximation order, one finds 

(29) KWV 1 =^t/ S r-V°, 

r 

so that 

(30) ( R+ l§i) ^ + ^ 2 )+ (^ + EE urs ) ^ = °- 

where the term U sr + U rs is recognized as the first order approximation of the "exchange 
energy". Since the wave function ip must describe the distribution of all the particles, 
a "multi-time functional" ip{t l , t a , t n ) can be defined, such that, as previously, 
ip(t, ■■■,t) = ip(t), so that the Schroedinger equation for the wave functional reads 



(31) 



(r s + v s (t s ) + ^(* x . = o. 



A reference frame can be found, where R s does not depend on time explicitly; here, the 
eigenvalue equation 

(32) f(Rs)=u Vl .. M J(h Vs ) 
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holds, u Ug being time-independent functions; the corresponding time-dependent functions 
are v Ul ... Un = u l / 1 ... l /„e _l ^-= Ust , which satisfy ([3Tj) at the 0"* order. The functions w are 
defined as functionals of t s , q s and the fields A, and their time dependence is given by 



(33) 



£■■•£« 



ULl>-\ ...LOn 



so that 
(34) 



= £-£■ 



ft is now possible to solve the system (|27|). so that, at the I s * order, the functional %jj 
reads 



(35) 



4' 1 = -J2[ R r + 



h_d_ 

~tdT r 



v r (t r )r 



where ^(t 1 ...t n ) = e~ 1 ^-^ Rrtyt ~ t ^ij/ ) {t) : it is easy to verify that -0° must be a linear 
combination of the eigenfunctions u vl ,,, Un , with time-dependent coefficients e~ l ^- Urt . 
One is therefore interested only in the I s * order time-independent matrix elements of the 
operator 



(36) 



JJaif — 



J t 1 =t 2 = ...=t 



which are found by considering the integral 



(37) 



h d 



dq 1 ... I dq n j dAvt^ [R r + ) V s (t S )V r (t r )v^ 



because of ([M]), the condition {v' s - uj s ) + (y' r - uo r ) + J2 m =£ s ,r(^'m ~ w ™) = must 
be fulfilled. After standard manipulation, one finds that the time-independent matrix 
elements are given by the operator 



(38) 



dfe^ tr - t) \y,(t),Vr(t R ]e- iR *t tr - t) , 



where 
(39) 



V a 



x(A*(x),J a (x)) + 



f dA-(x) 
V dx 



S s (x) +c.c. + 0(A 2 



the terms in A 2 must be neglected, since, at this order, no quantity has been developed 
up to higher powers of the field. This way, the quantities J s and S s do not depend on 
the fields any more, and commute with them. In particular, one finds for the current J s 
the formal solution 



(40) 



Js(y) = e <=h J s (y) e r ^fr- . 



so that the interaction operator reads 



(41) 



V* 



J a {x)A r {xY + [S a {x) 



dArjx) 
dx 
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where 
(42) 

A T (x) = ^ dy J d 3 x J r (y)D(x -y)+ (s r (y), ^g^^ ) = dy J d 3 y J^ f (y)D(x 

with J e ff defined as previously. 

The generalization of this formalism to a "many-component" field [4] (what we would 
call a vector field) can be accomplished via the substitution of the field A with the field 
Ai, i = 0, 1, 2, 3, so that, for example, the scalar product A* A is replaced by J^i e i A i A i, 
where eo = — 1> £ i,2.3 = 1, and so on. New commutation relations have to be introduced, 
such as 

(43) [At(x),A 3 (y)] = 2ihce i 6 ij D{x - y). 

In order to have a positive-definite energy density for the field Ai, the new field B, 
the so-called Stueckelberg field, has to be introduced; in fact, the energy density H = 
Y,i e i^( A i)^ where 



n{A) -^ \ \ dx^dx~ k 



( 44 ) nA) = i-\y:^-^-+fA*A 



contains a negative term when i = 0. The way followed by Stueckelberg in order to 
determine this term is the "Dirac-Fock-Podolski approximation" (^,v4) i/j — 0, which 
can be interpreted as a Gupta-Bleuer condition on the divergence of A: when a mass 
term is introduced, and when a vector field is taken into account, the approximation 
reads 

(45) - ^ ^-ip = (-divA*divA - l(B*divA + divA'B) - 1 2 B*b) tb 

oxq dxo V / 

that eliminates the negative term in (|44jl . After standard manipulation the energy density 
for the two fields, H'(A, B) reads 

8ttH' = (rotA*,rotA) + (gradA^ + ^, gr adA + §£) + 

( lA o - f£) ( lA o - H) + (l A *+3radB*) (lA + gradB) : 

if the new potential <f>i = Ai + ^il~ x §^- is introduced, the energy desity rewrites 

(46) U ' = h i( P *' P ) + (C*,G)+1 2 {$*,$) +f4>Uo] , 

where F = {F (n , F 02 F 03 } and G = {F 23 , F31F12}, F i2 being the field strength, Fij = 

d±j_ _ dfa = OAk dAj 
tl dxi t 3 dx 3 — fc * dx z t J dxj ■ 

If the same calculation as the case of the scalar field is followed, motion equations for 
spinors and bosons are obtained. E.C.G. Stueckelberg, in fact, wanted set up a unifying 
theory for scattering and decay processes, within the framework of boson "gauge" fields: 
he achieved this task by taking into account the then-known particles and interactions, by 
hypothesizing generalized- "charge" conservation, and by predicting, from his calculation, 
the existence of new particles and information about their masses [5 . Unfortunately, not 
all leptons had already been observed yet, and, consequently, the notion of leptonic and 
barionic number, as well as the distinction of weak and strong interactions, had not 
already been introduced at that time, but he laid the theoretical foundation of gauge 
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theories. According to the results of the experiments, he classified the known "spinor" 
particles according to their scattering and decay properties by attributing them electric 
and "heavy" charges, so that electrons, neutrinos, protons and neutrons are referred to 
as e(l, 0), n(0, 0), P(l, 1) and N(0, 1), respectively. As in modern gauge theories, inter- 
action between these charges are described by boson fields, which follow directly from 
the eigen-value equations for the generators, so that he predicts four such fields, that, 
according to the interaction they carry, are classified as e (1, 0), n(0, 0), P (1, 1) and N 
(0, 1), respectively. Therefore, the processes mediated by these fields are 



• 1) processes mediated by n(0,0): 

the only processes described by these fields are of the type 

(47) S^S' + n(0,0), 

where S can be referred to any kind of spinor. In Stueckelberg's interpretation, 
a better understanding of the proton-proton and neutron-neutron interactions in 
atomic nuclei could be achieved by means of the real field n(0, 0). 

• 2) processes mediated by the field e(l,0): 
these processes are /3-decays: 

(48) P(l,l) ~jV(0,l) + e(l,0), 
together with 

(49) e(l,0)~(-n(0,0)) + e(l,0) 

describe a nuclear decay, where the notion of anti-particle follows from the mathe- 
matics of the model. 

• 3) processes mediated by N (0, 1): 
the reactions 

JV(0,1) <- (-n(0,0)) + N(0,l), 



(50) P(l, 1) «-» e(l, 0) + N(0, 1) 

lead to estimate the mass of the particle N (0,1): since the proton is a stable 
particle, the mass of N (0, 1) must be greater than the difference of the masses of 
the proton and of the electron; furthermore, because of statistics, the mass of the 
particle N (0, 1) must be greater than the neutron mass, and it must be an instable 
particle, whose decay mode is 

(51) N(0,l)^P(l,l) + (-e(l,0)). 



• 4) processes mediated by P (1, 1): 



(52) P(l,l) «-» (-n(0,0))+P(l,l), 

so that the mass of the particle P(0, 1) must be greater than the proton mass. 

A modern approach to the Electroweak model via the Stueckelberg field is proposed in 
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4. The Stueckelberg field beyond Stueckelberg 

Application of this formalism was at first devoted to demonstrate its renormalizabil- 
ity. In this sense, Zimmermann[9] started to study Stueckelberg Lagrangian (|12p . and its 
invariance under the Pauli gauge transformations. At the end, the Stueckelberg massive 
Abelian model was proved to be renormalizable and unitary by Lowenstein and Schroer 
in 1972 [10J . We want to stress that this implies that the Stueckelberg model is the only 
way to give a mass to an Abelian boson, without a spontaneous symmetry breaking mech- 
anism. Therefore, there were several attempts to apply the theory to the non-Abelian 
case, in order to to furnish an alternative to the Higgs boson in the Standard Model. 
In 1988, Delbourgo, Twisk and Thompson^llJ first proved that the original Stueckelberg 
theory for neutral massive vector fields is invariant under nilpotent BRST transforma- 
tions, which ensures unitarity and renormalizability. Their work clearly illustrated that 
the key point, to avoid divergences, is the invariance under Pauli transformations. Then, 
they also analyzed the extension to non-Abelian fields. They noticed that renormaliz- 
ability and unitarity seem to be competing qualities of massive non Abelian theories, so 
they argued: 

"Finally, it must be admitted that the Higgs mechanism remains the most 
complete method for giving mass to the vector bosons" . 

But extension of the Standard Model, such to contain a Stueckelberg field, are again 
under investigationpj] . 

However, Stueckelberg theory for massive bosons found application also very far from its 
natural context. 

An example is given by the work of Ramond[13 , who applied the scheme to obtain 
the fully covariant and gauge invariant field theory for free open bosonic strings in 26 
dimensions. In fact, Stueckelberg fields naturally arise and are shown to be unrestricted 
for the most general gauge transformations. 
To quote his own words: 

"It should be clear that Stueckelberg field leads to much simpler looking ex- 



Moreover, also very recent attempts to introduce a massive Abelian field in Loop Quan- 
tum Gravity deal with the Stueckelberg field. 

Hence, HelesfaijM] stressed how, in such a context, the application of the Stueckelberg 
formalism is very useful since no second class constraint arise and the Hamiltonian is a 
linear combination of constraints (after quantization, the Proca field leads to a Hamilto- 
nian that is quadratic in the Lagrange multipliers). In fact, the Hamiltonian reads 



H = ±tr(2[K a ; K b ] F ab )[E«;E b ] + jf q (E a E b + B a B b ) + ^ + ^q ab (A a + 
da<f){A b + d b 4>) 

K a = F 3 ab E b + e abc E b B c + (A a + o )tt 



pressions" . 



(53) 




being 
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G = D a E a - 7T 
Gb = D a E^ 

The quantization is performed on the Hilbert space 

(54) H = L 2 {A S u(2),dnsu(2)) ® L 2 {A u{1) ,d^ u{1) ) <g L 2 {U u{1) ,d^ u(1) ) 
for which a basis is given by the generalized spin network functions 

( 55 ) \S > 7J V, J >= \ T (A) > 7tl ?®\F(A) > iX ®\D(U) > lM ■ 

In this context, the mass m is a coupling constant and is very similar to the Immirzi 
parameter (in the quantum regime, it enters the Hamiltonian in a non-trivial way). 

5. Brief concluding remarks 

Among the brilliant results accomplished by Stueckelberg, the formulation of a divergence- 
free model for massive vector fields has been one of the most prolific ideas in modern 
Physics. In fact, despite the Proca formulation, his intuition of the need to maintain a 
gauge invariance in the theory has been the key to the latcr-rccognizcd rcnormalizability. 
Moreover, the modernity of his approach relies on the preference of a gauge symme- 
try rather than phenomenological speculations, such as the Yukawa formulation [2], and 
renders the Stueckelberg field a suitable tool also in current achievements of theoretical 
Physics, i.e. String theory and LQG. 
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